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ABSTRACT 

We model the atmospheres and spectra of Earth-like planets orbiting the entire grid of M dwarfs for 
active and inactive stellar models with T e s = 2300K to T e s = 3800K and for six observed MUSCLES 
M dwarfs with UV radiation data. We set the Earth-like planets at the 1AU equivalent distance and 
show spectra from the VIS to IR (0.4/rm - 20/xm) to compare detectability of features in different 
wavelength ranges with JWST and other future ground- and spaced-based missions to characterize 
exo-Earths. We focus on the effect of UV activity levels on detectable atmospheric features that 
indicate habitability on Earth, namely: H 2 0, O3, CH4, N 2 0 and CH 3 C1. 

To observe signatures of life - O2/O3 in combination with reducing species like CH4, we find that 
early and active M dwarfs are the best targets of the M star grid for future telescopes. The 0 2 
spectral feature at 0.76^m is increasingly difficult to detect in reflected light of later M dwarfs due 
to low stellar flux in that wavelength region. N 2 0, another biosignature detectable in the IR, builds 
up to observable concentrations in our planetary models around M dwarfs with low UV flux. CH 3 C1 
could become detectable, depending on the depth of the overlapping N 2 0 feature. 

We present a spectral database of Earth-like planets around cool stars for directly imaged planets as 
a framework for interpreting future lightcurves, direct imaging, and secondary eclipse measurements 
of the atmospheres of terrestrial planets in the HZ to design and assess future telescope capabilities. 

Subject headings: astrobiology, planets: atmospheres, planets: terrestrial planets, stars: low mass 


1. INTRODUCTION 

About 2000 extrasolar planets have been found to date 
with thousands more awaiting confirmation from space 
and ground-based searches. Several of these planets 
have been found in or near the circumstellar Habitable 


Kaltenegger et al. 2013; Bataiha et al. 

2013; Borucki 

et al. 2011: IKaltenegger & Sasselov 201 
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models. Future mission concepts to characterize Earth¬ 
like planets are designed to take spectra of extrasolar 


mospheric signatures (e.g. Beichman et al. 1999 

2006 

Cash 

2006 

Traub et al. 

2006 Kaltenegger et al. 

2006) 


nearby Super-Earth and Earth-like planets using emer¬ 
gent visible and infrared spectra. For transiting terres- 
trial planets, the James Web Space Telescope (JWST, 


see Gardner et al. 12 006] |Deming et al.|2 009| [Kaltenegger 

&; Traub|2009D as well asTuture ground ana space based 


telesc opes (ISnellen et al.||2013[ |Rodler & Lopez-Morales 
20141 will search tor biosignatures in a rocky planet’s at- 
mosphere. NASA’s explorer mission, TESS, is designed 
to search the whole sky for potentially habitable p lanets 
around the closest and brightest stars to Earth (et al. 
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20141 for eventual follow-up with JWST and other large 


ground-based observatories such as the E-ELT or GMT. 

In our solar neighborhood, 75% of stars are M dwarfs. 
The abundance of M dwarfs as well as the contrast ra¬ 
tio and transit probability favor the detection of plan¬ 
ets in the Habitable Zone of M dwarfs. Therefore, it is 
likely that the first habitable planet suitable for follow¬ 
up _obseryations_will_beJound_orbiting a nearby M dwarf 
(Dressing & Charbonneau 2013). The M spectral class is 
very diverse, spanning nearly three orders of magnitude 
in luminosity and an order of magnitude in mass. 

The UV environment of a host star dominates the pho¬ 
tochemistry and therefore the resulting atmospheric con¬ 
stituents including biosignatures for terrestrial planets 


(see e.g. Domagal-Goldman et al. 2014 

Tian et al. 2014 

Grenfell et al. 20141 IRughelmer et al. 

2013; Hu et al. 

2012| Segura et aTpOOSj). To date, few 

observations ex- 


spectra of very active flare stars, such as AD Leo, were 
available from the IUE satellite, primarily during flares 
and for a few quiescent phases. The MUSCLES program 
observed chromosph eric emission from six weakly active 
M dwarfs with HST (France et al.|2013|. Note that Ly-a 
is by far the strongest lineTn the UV for M dwarfs. Since 
the core of the intrinsic stellar Ly-a emission is absorbed 
by neutral hydrogen in the interstellar medium, one must 
reconstruct the line to compens ate for this absorption to 
get accurate flux levels (see e.g. iWood et al. 2005 Linsky 

Several groups have explored the effect of a different 
stellar spectral type on the atmospheric composition of 
Earth-like planets by primarily considering one star, AD 
Leo, as a template for the diverse range of M dwarfs, 
and by using the extreme limit of inactivity, photosphere 
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only, PHOENIX models ([Grenfell et al. [20141 Kitzmann 
et al.||2011a|bl |Wordsworth et al | |2011| |Segura et al. 
201 ).V;. In t his paper we expand on this work by es¬ 
tablishing planetary atmosphere models for the full M 
dwarf main sequence, using a stellar temperature grid 
from 3800K to 2400K including recent HST observations 
of 6 M dwarfs to explore the effect of different spec¬ 
tral energy distributions on terrestrial atmosphere mod¬ 
els and on detectable atmospheric signatures, including 
biosignatures. Atmospheric biosignatures are remotely 
detectable chemical species in the atmosphere that are 
byproducts of life processes. Recent research shows pos¬ 
sible fals e posit i ves may occ ur under c ertain geolog ical 
settings ( Domagal-Goldman et al . 1 2014 |Tian et al7|2014| 


Wordsworth & Pierrehum bert|| 2 (TM ). 


We create a grid of IVl dwarf input spectra, that can 
be used to probe the entire range of UV flux levels. To 
showcase the whole range, we simulate model planets 
around host stars at the two extreme limits of activity: 
active and inactive stellar models. We compare these 
models to th e six observed M dwarfs with recent UV 
observations (France et al.||2013). 


We explore the influence oFstellar UV flux on the at¬ 
mospheric structure, chemical abundance, and spectral 
features for Earth-like planets including the observability 
of biosignatures in the VIS to IR. We focus our analysis 
on spectral biosignatures for a temperate rocky planet 


fluence of stellar types on the abundance of atmospheric 
chemical species. In §4, we examine the remote observ¬ 
ability of such spectral features, and in §5 and § 6 , we 
summarize the results and discussing their implications. 

2. MODEL DESCRIPTION 
2 .1. Stellar M Dwarf Spectral Grid, Model 


represent the most active M dwarf measurements and 
“inactive” semi-empirical models without chromospheres 
are constructed to represent the lowest theoretical UV 
flux field. We compare these limiting-case models to six 
well-observed M dwarfs which show significant chromo¬ 
spheric flux (see Fig. [lj despite being traditionally clas¬ 
sified as qu iescent stars due to the presence of H Q in 
absorption ( France et al.|[2013| ). A few M dwarfs ob¬ 
served with GALEX in the NUV (1750A- 2750A) show 
near photospheric continuum level fluxes (as computed 
with PHOENIX models) and can be classified as inactive 
even though the exact UV flux level has not been mea¬ 
sured yet. For example, G1 445 and G1 682 have roughly 
90% of the observed near ultraviolet (NUV) flux d ensity 
is predicted by PHO ENIX model photospheres (|Shkol-| 
nik & Barman 2014). More observations are needed to 
determine the ultimate lower limit of UV flux emitted by 
old and/or late M dwarfs. To compare our calculations 
for inactive star models to published work, we use theo¬ 
retical photos phere PHOENIX models as a lower bound 
( Allard 2~0l4|). 
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M dwarfs span nearly three orders of magnitude in lu¬ 
minosity and remain active for much longer timescales 
than earlier stellar types. As main sequence stars age, 
their UV flux levels decrease even as their bolometric lu¬ 
minosity increases. Understanding the UV flux incident 
on a planetary atmosphere is critical to understanding 
and interpreting future observations of atmospheric con¬ 
stituents, including biosignatures. 

Current stellar models are unable to model the UV 
region from M dwarfs self-consistently due to three main 
reasons. First, complex magnetic fields responsible for 
heating the chromosphere are thought to drive much of 
the UV activity and have thus far been ignored in stellar 
models. Second, the models are missing opacities in the 
UV, and third, the models are semi-empirical with no 
energy conservation to balance magnetic heating with 
radiative losses. These problems are being addressed by 
current work and models will be available in the future 
to test against M dwarf UV observations 

We generate stellar input models for the entire M dwarf 
spectral class (M0 to M9) to explore the boundaries of 
the UV environment of an exoplanet orbiting a main se¬ 
quence star. We create two sets of models based on the 
extreme limits of stellar activity. For the purposes of 
this paper, “active” stellar models are constructed to 


Figure 1. Stellar input spectra at the top-of-atmosphere (TOA) 
of Earth-like planet at the 1AU equivalent of an M0 “active” star 
with UV flux scaled from AD Leo (black), “inactive” PHOENIX 
model of M0, and HST UV observations plus photosphere model 
for GJ436 (green). 


The chromospheric far-ultraviolet (FUV) flux from M, 
L and T dwarfs (i.e., very low mass stars and brown 
dwarfs), while crucial for determining the potential hab¬ 
itability of any planets around them, is very poorly char¬ 
acterized through direct observations. Given their com¬ 
mon chromospheric origin, H a , Ly-a, Ca II H and K, 
and Mg II H and K emission lines have all been used as 
proxies for FUV activity. However, many of these lines 
are inaccessible to M dwarfs because they are at shorter 
wavelengths where these stars are intrinsically less lu¬ 
minous. For M dwarfs, H a emission in particular has 
been studied extensively and to date has the most ro¬ 
bust dataset for each M dwarf spectral ty pe co m pared 


with the other line s considered above (West et al.||2004 
2011|) . In addition, | Jones & West| ([2014 ) shows H Q fluxes 
correlate with NUV Huxes from GALEX. Therefore, we 
use H a to estimate the FUV emission for our active stel¬ 
lar grid by scaling the known FUV and H a emission from 
the active M3.5 star AD Leo (see Eq. [l] and Fig# Very 
active, early M dwarfs (M0-M5) are known to saturate in 
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the H„ emission around log(Lf/„ /Lh n i) = —3.75 (Hawley 

et al. 1996 

Reiners & Basri 2008 

West et al. 2004). 
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Figure 2. Scaling us ed for active stellar models (black line) over 
plotted on|West et all] p004j measurements of log(mean Ljj a /Lf, 0 i) 
(black asterisks). 
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West et al. (2004) tabulates the log(mean Ln^/Lboi) 
versus spectral type from M0-L0 from observations of 
1910 active M dwarfs. We parameterize the data from 
West et al. (2004|) to derive a relationship of UV to mean 


H a emission for each spectral sub type in the M dwarf 
class (See Fig. [2] and Eq[l|. 


log (L H ^/L boU ) = -3.62 for M0 - M4 (1) 

log(LH a */Lb 0 u) = — 0.32 x (M) — 2.07, 

where M = 5, 6 , 7 for M5, M 6 , M7 
\og{L Ha jL boU ) = - 4.25 for M 8 - M9 


We use the H a /UV scaling of AD Leo, which has both a 
well characterized UV spectrum and H a emission mea¬ 
surement as a calibration data point to set the scale. 
AD Leo, an M3.5eV star, has log (LH a /L bo i) = -3.51 
(Walkowicz &; Hawley 2009) an d has been well observed 
in the UV (s ee [begura et al.| (|20 05|) for AD Leo IUE 
spectrum and Linsky et al. (|2013 ) for reconstructed Ly- 
a flux). The M dwarf active spectrum is then given by: 


FfUV * = Fpuv AD Leo X ( 


Teff 


7 efF AD Leo 

. log {LH a */Lbou) 

log {Ln a ADLeo/LbolADLeo) 


( 2 ) 


Our stellar temperature grid (Fig. [3]) covers the full M 
dwarf spectral range with the stellar parameters for ac¬ 
tive and inactive models given in Table [II For the ac¬ 
tive star models, we scale the UV flux in the wavelength 
range 1000 - 3000A from AD Leo’s observed UV spec¬ 
trum to the other spectral types according to Eq. [2] For 
each active model star on our grid we concatenated a 
solar metallicity, unreddened synthetic PHOENIX s pec- 
trum, which only considers photospheric emission (Al¬ 


lard 20141, from 3000 - 45,450A to the scaled observa¬ 
tions from the International Ultraviolet Explorer (IUE) 
archive for AD Leo from 1000 - 3000A combined with 


the reconstructed Ly-ct from Linsky et al. (2013) for the 


region 1210 - 1222AQ Note that young active stars are 
highly variable, especially in the UV. AD Leo is the most 
active and well characterized M dwarf. Our active mod¬ 
els therefore represent the extreme high end of activ¬ 
ity for each stellar type from an active flare star. For 
comparison, the inactive semi-empirical models include 
no chromospheric emission. The continuum originating 
from the photosphere is taken from a PHOENIX model 
with the same stellar parameters as the corresponding 
active star from 1000 - 45,450A. Thus, these models pro¬ 
vide only a lower limit to the stellar UV flux. More ob¬ 
servations are needed to determine the true lower limit 
of inactive M dwarfs in the UV. 


Table 1 

Stellar properties for active and inactive models. 


Star 

Us (U) 

Mass ( M(7 )) 

Radius (ifo) 

M0 

3800 

air 

0.62 

Ml 

3600 

0.49 

0.49 

M2 

3400 

0.44 

0.44 

M3 

3200 

0.36 

0.39 

M4 

3100 

0.2 

0.36 

M5 

2800 

0.14 

0.2 

M6 

2600 

0.10 

0.15 

M7 

2500 

0.09 

0.12 

M8 

2400 

0.08 

0.11 

M9 

2300 

0.075 

0.08 


The models represent the extreme limits of the UV ra¬ 
diation environment for an exoplanet using a flare star 
for the upper bound and a photosphere-only model for 
the lower bound in the UV. Fig. [I] shows the UV flux 
region between those two limits spans over 10 orders of 
magnitude in the FUV. We use the observations of six M 
dwarfs from the MUSCLES program to probe the region 
in between our active and inactive models. The MUS¬ 
CLES stars have been traditionally classified as “quies¬ 
cent” because they do not present H a in emission. 

The six M dwarfs were observed by two UV spectro - 
graphs (COS and STIS) on HST ( [France et al. 2013). 
We joined these HST measurement^ with the star’s cor¬ 
responding PHOENIX models at 2800A after adjusting 
the HST flux levels to the level a planet would receive 
at the 1AU equivalent in the Habitable Zone. The input 
parameters for the PHOENIX models are the observed 
star’s T e ff, [Fe/H], log g , and the rotation velocity, v sin 
i, as summarized in Tableland described for each star 
based on the most recent observations in the Appendix. 
When observations of v sin i were not available we used 
instead the peak o f the observed distribution for 56 M 
dwarfs (Jenkins et al.||2009). 


The Muscles database interpolates the region be¬ 
tween 1760 - 2100 A because the noise in the observa¬ 
tions was large in that region since those wavelengths 
were covered by the lowest sensitivity part of the STIS 
G230L bandpass. In order to not bias the results too 
higher or too low, the MUSCLES team used the average 
flux in the 2100 - 2200 A region to interpolate the 1760 
- 2100 A region. Follow up observations to the MUS¬ 
CLES program with a new HST Treasury program aim 
to probe the 1760 - 2100 A region of the spectrum with 
higher sensitivity. We used the MUSCLES database and 

7 http://archive.stsci.edu/iue/ 

8 http://cos.colorado.edu/ kevinf/muscles.html 
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Table 2 

Stellar properties for MUSCLES model stellar spectra. 


MUSCLES 

Stars 

Teff (K) 

Radius (Rq) 

[Fe/Hj 

log 9 

Age (Gyr) 

v sin i 
(km s -1 ) 

GJ 832 

3620 

0.48 

-0.12 

4.70 

- 

3 * 

GJ 667C 

3350 

0.348 

-0.55 

5.00 

> 2 

3 * 

GJ 1214 

3250 

0.211 

+0.05 

4.99 

6+3 

< 1 

GJ 436 

3416 

0.455 

+0.04 

4.83 

6.5-9.9 

< 1 

GJ 581 

3498 

0.299 

-0.10 

4.96 

7-11 

2.1 

GJ 876 

Note: for GJ 832 

3129 

and GJ 667 

0.3761 

CJ there is no v sin i 

+0.19 

measureme 

4.89 

mt. We ; 

0.1-5 

assumed a v sin i ■■ 

1.38 

= 3.U corresponding 


to the peak of the distribution {Jenkins et al.|2009) when generating a PHOENIX model. 
See text for references. 


interpolation as given. Because this region of the UV is 
important for O3 production, uncertainty of this inter¬ 
polation could influence our atmospheric results. 

Input active stellar spectra and MUSCLES spectra are 
shown in Fig. [3] (inactive models are not shown). 

2.2. Planetary Atmosphere Model 
We use EXO-P ( Ealtenegger &; Sasselov| [2010), a cou¬ 


pled ID radiative-convective atmosphere code developed 


1 Easting & Ackerman 

1986 

Pavlov et al.; 2000-; Haqq- 

Misra et aT.T2l)U8 j. ID p 

tiotochemistry (Pavlov &; Pasting 

2002 Segura et al. 2005 

:2007), and ID radiative transfer 

model (Traub & Stier|1976 

Kaltenegger & Traub 2009) 


planet orbiting host M dwarfs in the Habitable Zone. 

EXO-P is a model that simulates both the effects of 
stellar radiation on a planetary environment and the 
planet’s outgoing spectrum. We model an altitude range 
that extends upwards to 60km with 100 height layers. 
We use a geometrical model in which the average ID 
global atmospheric model profile is generated using a 
plane-parallel atmosphere, treating the planet as a Lam¬ 
bertian sphere, and setting the stellar zenith angle to 
60 degrees to represent the average in coming stellar flux 


on the da yside of the planet (see also Schindler & East¬ 


ing 2000). The temperature in each layer is calculated 


from the difference between the incoming and outgoing 
flux and the heat capacity of the atmosphere in each 
layer. If the lapse rate of a given layer is larger than 
the adiabatic lapse rate, it is adjusted to the adiabatic 
rate until the atmosphere reach es equil ibrium. We use a 


two-stream approximation (see Toon et aT7|| 1989), which 


includes multiple scattering by atmospheric gases, in the 
visible/near IR to calculate the shortwave fluxes. Four- 
term, correlated-k coefficients pa rameter i ze the ab sorp- 
tion by O3, H 2 0, 0 2 , and CH4 {Pavlov et al.||2000|). In 
the thermal IR region, a rapid radiative transfer model 
(RRTM) calculates the longwave fluxes. Clouds are not 
explicitly calculated. The effects of clouds on the temper¬ 
ature vs. pressure profile are included by adjusting the 
surface albedo of the Earth -Sun system to have a surface 
temperature of 288K (se e Easting et al. 1984 Pavlov 
et al.||2000| ISegura et a.l.||2()03[ |2()( io[). The photochem- 
lstry code, originally developed by|Rasting et al.| {1985 1 
solves for 55 chemical species lin ked by 221) reaction s us¬ 
ing a reverse-Euler method (see Segura et al.||2010 and 
references therein). The photochemical model is station¬ 
ary, and convergence is achieved when the following cri¬ 
teria are fulfilled: the production and loss rates of chem¬ 
ical species are balanced which results in a steady state 
for the chemical concentrations, and the initial boundary 


conditions, such a surface mixing ratios or surface fluxes, 
are met. 

The radiative transfer model used to compute plane¬ 
tary spectra is based on a model originally develo ped fo r 
trace gas retr ieval in Earth’s atmospheric spectra (Traub 
fe Stier|1976 ) and further devel oped for exoplanet trans¬ 
mission and emergent spectra (Ealtenegger et al.||2007 


Ealtenegger & Traub 2 0091 |EalteneggerI2 010; Ealteneg- 

ger & Sasselov||2010| {Raltencgger et ah||2013 k In this 

paper, we model Earth's reflected and thermal emission 
spectra using 21 of the most spectroscopically significant 
molecules (H 2 0, O3, 0 2 , CH4, C0 2 , OH, CH3CI. N0 2 , 
N 2 0, HN0 3 , CO, H 2 S, S0 2 , H 2 0 2 , NO, CIO, HOC1, 
H0 2 , H 2 CO, N 2 0 5 , and HC1). We use a Lambert sphere 
as an approximation for the disk integrated planet in our 
model. The surface of our model planet corresponds to 
Earth’s current surface of 70% ocean, 2% coast, and 28% 
land. The land surface consists of 30% grass, 30% trees, 
9% granite, 9% basalt, 15% snow, and 7% sand. Surface 
reflectivities are taken from the USGS Di gita l Spectral 
Lib ranFl and the ASTER Spectral Library^] (following 
Ealtenegger et al. 2007). 

Clouds have a strong impact on the detectability of at¬ 
mospheric species. For the spectra shown in Figs. m 
we assume a 60% global cloud cover with cloud layers 
analogous to Earth (40% water clouds at 1km, 40% wa¬ 
ter clouds at 6km and 20% ice clouds at 12km follow¬ 


ing Ealtenegger et al. 2007). In the VIS to NIR, clouds 
increase the reflectivity of an Earth-like planet substan¬ 
tially and therefore increase the equivalent widths of all 
observable features, even though clouds block access to 
some of the lower atmosphere. In the IR, clouds slightly 
decrease the overall emitted flux of an Earth-like planet 
because they radiate at lower temperatures and therefore 
decrease the equivalent widths of all observable absorp¬ 
tion features, even though they can increase the relative 
depth of a spectral feature due to lowering the continuum 
temperature of the planet. For a comparison of scenarios 
wit h Earth-analogue clouds to those of clear sky spectra 


see 


Rugheimer et al. (2013). 


Using 34 layers, we calculate the spectrum at high 
spectral resolution with several points per line width. 
The line shapes and widths are computed using Doppler 
and pressure broadening on a line-by-line basis for 
each layer in the model atmosphere. The overall 
high-resolution spectrum is calculated with 0.1 cm -1 
wavenumber steps. The figures are shown smoothed to 
a resolving power of 150 in the IR and 800 in the VIS 
using a triangular smoothing kernel. The spectra may 

9 http://speclab.cr.usgs.gov/spectral-lib.html 

10 http://speclib.jpl.nasa.gov 
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Figure 3. Stellar input spectra (top) from 1000 to 20000 A for the M0-M9 active grid stars with UV scaled by H a and AD Leo UV flux 
(left) and MUSCLES stars with HST UV observations (right). UV input stellar fluxes in log scale (bottom). 

be binned further for comparison with proposed future 
spectroscopy missions designs to characterize Earth-like 
planets. We previously validated EXO-P from the VIS 
to the IR u sing data from ground and space (IKaltenegger 

gne 


= 9.0 x 10 8 , and CCH 3 CI = 5.0 x 10 10 (see Rugheimer 


et al.||2007[ IKaltenegger & Traub 2009 Rugheimer et al. 

20T3f. 


2.3. Simulation Set-Up 

To examine the effects of the varying UV flux of M 
dwarfs on an Earth-like atmosphere and its observable 
spectral features, we use the temperature grid of stel¬ 
lar models ranging from M9 to M0 ( T e g = 3800K to 
2300K) for active and inactive models along with the 
six M dwarfs with well characterized UV fluxes from 
HST (see §2.2). We simulated an Earth-like planet with 
the same mass as the Earth at the 1AU equivalent or¬ 
bital distance, defined by the wavelength integrated stel¬ 
lar flux received on top of the planet’s atmosphere be¬ 
ing equivalent to 1AU in our solar system calculated by 
1 AU eq =AUx J{R/Rq) 2 X (T eff /T eff 0 ) 4 . 

The biogenic (produced by living organisms) fluxes 
were held fixed in the models in accordance with the 
fluxes that reproduce the modern mixing ratios in the 
Earth-Sun case, except for the cooler M dwarfs where 
CH 4 and N 2 O were given a fixed mixing r atio of 1 .0 x 
1CU 3 and 1.5 x 10 —2 , respectively (following Segura et al. 


2003 2005). The surface fluxes for the long-lived gases 
H 2 , CH 4 , N 2 0, CO and CH 3 C1 were calculated such 
that the Earth around the Sun yields a T sur f = 288K 
for surface mixing ratios: cH 2 = 5.5 x 10 -7 , CCH 4 = 
1.6 x 10 -6 , cC0 2 = 3.5 x 10" 4 , cN 2 0 = 3.0 x 10" 7 , cCO 


et al. 2013). The corresponding input surface fluxes to 
the atmosphere are —1.9 x 10 12 g H 2 yr _1 , 5.3 x 10 14 g 
CH 4 yr-\ 7.9 x 10 12 g N 2 0 yr" 1 , 1.8 x 10 15 g CO yr" 1 , 
and 4.3 x 10 12 g CH 3 CI yr _1 . The N 2 mixing ratio is 
set to be a “fill gas” such that the total surface pressure 
is 1 bar. These boundary conditions were used for MO¬ 
MS active grid stars, M0-M3 inactive grid stars and all 
MUSCLES stars. 

For M6-M9 active grid stars and for M4-M9 inactive 
grid stars, the boundary condition for CH 4 is changed 
to a fixed mixing ratio once the UV environment of the 
host star drops below a certain level because, CO, H 2 , 
and CH 4 do not converge according to the criteria de¬ 
scribed i n §2.2 assu ming a modern Earth biological flux 


(see also Segura et al. 


20051. In reality, the thermody¬ 


namical cost for microbes producing CH 4 would become 
unprofitable as temperatures would initially rise due the 
increased greenhouse effect, causing the microbes to op¬ 
erate less efficiently. This type of biological feedback 
is not included in the models. CH 4 is produced biot- 
ically by methanogens and other organisms and abioti- 
cally through hydrothermal vent systems. In the modern 
atmosphere there is a significant anthropogenic source of 
CH 4 from natural gas, livestock, and rice paddies. Up¬ 
per limits of abiotic fluxes of methane can be estimated 
f or ter restrial planets following Guzman-Marmolejo et al. 
(2014|. An estimate of Earth bio tic m ethane is abou t 30x 
the abiotic flux (see discussion in Segura et al.|2005|), and 
thus a range of methane fluxes may be maintained on a 
terrestrial planet via methanogenesis in addition to abi- 
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otic methane production. 

Therefore, we set the mixing ratio of CH 4 to 1 x 10 -3 
corresponding to a value of the last stable run with an 
Earth-like biological CH 4 flux for an active M5 star for 
planets around host stars that show this “run-away” be¬ 
havior. Given these boundary conditions, the methane 
flux necessary to sustain a 1 x 10 -3 mixing ratio of 
methane is 5.65 x 10 14 g yr _1 (equal to our present Earth 
methane flux) for the planet around the hottest M dwarf 
(T e g = 3800 K) either active or inactive, 1.16 x 10 14 g 
yr _1 (20.5% of the Earth value) for the planet around 
the coolest (T e g = 2300 K) active star, and 5.87 x 10 12 g 
yr _1 (1% of the Earth value) for the coolest inactive M 
dwarf of our sample. 

For the less active stars, we set the deposition velocity 
(rate at which they are deposited to the surface) for H 2 
and CO to 2.4 x 10 -4 cm s _1 and 1.2 x 10 -4 cm s _1 , 


respectively (see also Segura et al. 20051 

Kharecha et al. 

2005) 

Rauer et al. ( 

2011 

) used a deposil 

cion velocity tor 


after calibrating the model to reproduce the H 2 atmo¬ 
spheric abundance measured for present Earth. In our 
case, the deposition values for H 2 and CO correspond to 
maximum air-sea transfe r rates estimated using the “pis¬ 


ton v elocity” approach (Broecker 1982 Khareclra et al. 


20051. Since a larger H 2 deposition velocity implies con- 


sumption by ba cteria , we use the abiological value in 
Kharecha et al. (20051. 

We followed a similar procedure for constraining the 
N 2 O concentrations for M6-M9 inactive stars. We set 
the mixing ratio of N 2 O to 1.5 x 10 -2 , corresponding 
to values of the last stable run with Earth-like biologi¬ 
cal fluxes, to prevent unphysical buildup of N 2 O. N 2 O is 
emitted primarily by (de) nitrifying bacteria with anthro¬ 
pogenic sources from fertilizers in agriculture, biomass 
burning, industry, and livestock and is a relatively mi¬ 
nor constituent of the modern atmosphere at around 320 
parts per billion (ppb), com pared to the^ pre-industrial 
concentration of 270 ppb (Forster et al. 2007). 

We used modern Earth fluxes for ail of the MUSCLES 
stars without any “run-away” buildup of reduced gases. 
The CH 4 and N 2 O fluxes are given in Table [3] 

All of our simulations used a fixed mixing ratio of 
355ppm for CO 2 , 21 % O 2 , and a fixed upper boundary of 
10 -4 bar (~ 60km). In the ID climate model, a surface 
albedo of 0.2 is fixed for all simulations, corresponding to 
the surface albedo that reproduces Earth’s average tem¬ 
perature of 288K for the Earth/Sun case. The planetary 
Bond albedo (surface + atmosphere) is calculated by the 
ID climate code.cm -2 


3. ATMOSPHERIC MODEL RESULTS 

The amount of UV radiation emitted from the host 
star influences the abundances of major chemical atmo¬ 
spheric constituents including H 2 O, CH 4 , and O 3 and, 
as a result, modified the temperature-pressure profile of 
a planet. The UV fluxes incident at the top of the at¬ 
mosphere of an Earth-like planet in the HZ are given in 
Table [4] and shown in Fig. [3] An M0 active star model 
has 6.7 times more total UV flux than an M0 inactive 
star model. An M9 active star model has 4600 times 
more total UV flux than an M9 inactive star model. The 
greatest differences are in the far UV (FUV 1000 - 2000 


A) where an active star model has 6.8 x 10 4 and 1.3 x 10 11 
more FUV flux than an inactive M0 and M9 model, re¬ 
spectively. The MUSCLES stars’ UV environments fall 
between the active and inactive star models, consistent 
with their classification as being weakly active. 

The temperature vs altitude profile and the H 2 0, 0 3 , 
CH 4 and N 2 O mixing ratio profiles for all of the simula¬ 
tions are shown in Fig. El with each row corresponding 
to active, inactive, and MUSCLES star models, respec¬ 
tively. CH 3 CI profiles are not shown here, but follow 
the same trends as CH 4 . Since both 0 2 and C0 2 are 
well mixed in the atmosphere, their vertical mixing ratio 
profiles of 0.21 and 355ppm, respectively, are not shown. 

In the first column of Fig. El we show the changes 
in the temperature/altitude profile for Earth-like atmo¬ 
sphere models around M dwarfs for active (top row), in¬ 
active (middle row) and the six observed MUSCLES stars 
(bottom row). All temperature inversions are weaker 
than for the modern Earth because M dwarfs emit low 
UV flux in the 2000 - 3000 wavelength region, thereb y 


pro d ucing n ear isothermal stratospheres (see also Segura 
et al. 20051. Fig. El shows that temperature inversions 


are weaker for the higher UV environment stars. This is 
counter-intuitive since for the modern Earth, O 3 absorp¬ 
tion of UV radiation causes stratospheric heating and 
an inversion. However in these lower UV environments, 
additional heating is provided by stratospheric CH 4 and 
H 2 0. 

H 2 0 concentrations are lower in the stratosphere for 
planets around M dwarf models with higher UV fluxes, 
although higher O 3 concentrations will act as a shield, 


partially offsetting the effect of higher UV photon fluxes 
(Fig. El 2 nd column). H 2 0 can also be formed in the 
stratosphere from CH 4 and OH (CH 4 + OH -A CH 3 + 
H 2 0) and by increased upwards vertical transport in the 
nearly isothermal stratos pheres o f Earth -like planets or¬ 
biting M dwarfs (see also JSegura et a l. 2005). O 3 shields 
H 2 0 in the troposphere from U V environments. While 
photochemically inert in the troposphere, H 2 0 can be re¬ 
moved by photolysis at wavelengths shortward of 2000A 
in the stratosphere (see columns 2 & 3 in Table [4]) or 
by reactions with excited oxygen, 0( 4 D), to produce OH 
radicals. 

In an atmosphere containing 0 2 , O 3 concentrations 
are determined primarily by the absorption of UV light 


shortward of 2400A in the stratosphere, and we see a cor¬ 
responding decrease in O 3 concentrations correlated with 
decreased FUV radiation for cooler M dwarfs as well as 
inactive versus active stars (Fig. [4j 3rd column). 

CH 4 mixing ratio profiles are shown in column 4 of Fig. 
[4j As mentioned previously, we set the mixing ratio of 
methane to be 1000 ppm for the later M dwarf spectral 
types as described in §2.3 and Table [3j For Earth-like 
CH 4 fluxes, CH 4 concentrations decrease with higher UV 
environments due to reactions with OH to form H 2 0 and 
by photolysis by A < 1500 A (Fig. [ZJ 4th column). 

N 2 0 mixing ratios are larger for cooler M dwarfs be¬ 
cause cooler stars have smaller UV fluxes at A < 2200 
Aand thus lower photolysis rates. Around later M 
dwarfs, we see mixing ratios of N 2 0 that increase unre¬ 
alistically for the inactive stellar models M6-M9, similar 
to CH 4 (Fig. [4j 5th column). Since such a “run-away” 
effect should not be sustainable by biology, we cap the 
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Table 3 

CH4 and N2O fluxes and mixing ratios and O3 column depths. 



Surface 
Mixing 
Ratio CH4 

Flux (g/yr) 

ch 4 

% Earth 
CH 4 flux 

Surface 
Mixing 
Ratio N2O 

Flux (g/yr) 

n 2 o 

% Earth 
N2O flux 

O3 Column 
Depth 
(cm -2 ) 

Active 

MO A 

330ppm 

5.65 x 10 14 

100% 

0.70ppm 

7.99 x 10 12 

100% 

4.1 x 10 18 

Ml A 

370ppm 

5.65 x 10 14 

100% 

0.70ppm 

7.99 x 10 12 

100% 

4.0 x 10 18 

M2 A 

450ppm 

5.65 x 10 14 

100% 

0.76ppm 

7.99 x 10 12 

100% 

3.8 x 10 18 

M3 A 

570ppm 

5.65 x 10 14 

100% 

0.82ppm 

7.99 x 10 12 

100% 

3.6 x 10 18 

M4 A 

910ppm 

5.65 x 10 14 

100% 

0.93ppm 

7.99 x 10 12 

100% 

3.5 x 10 18 

M5 A 

1000 ppm 

5.65 x 10 14 

100% 

0.98 ppm 

7.99 x 10 12 

100% 

2.6 x 10 18 

M6 A 

lOOOppm 

2.90 x 10 14 

51.8% 

1.7ppm 

7.99 x 10 12 

100% 

3.1 x 10 18 

M7 A 

lOOOppm 

1.28 x 10 14 

22.7% 

3.0ppm 

7.99 x 10 12 

100% 

2.4 x 10 18 

M8 A 

lOOOppm 

1.36 x 10 14 

24.0% 

3.1 ppm 

7.99 x 10 12 

100% 

2.5 x 10 18 

M9 A 

lOOOppm 

1.20 x 10 14 

21.2% 

3.5ppm 

7.99 x 10 12 

100% 

2.3 x 10 18 

Inactive 

MO I 

490 ppm 

5.65 x 10 14 

100% 

34 ppm 

7.99 x 10 12 

100% 

1.7 x 10 18 

Ml 1 

580 ppm 

5.65 x 10 14 

100% 

64 ppm 

7.99 x 10 12 

100% 

1.5 x 10 18 

M2 I 

650 ppm 

5.65 x 10 14 

100% 

120 ppm 

7.99 x 10 12 

100% 

1.4 x 10 18 

M3 I 

1000 ppm 

5.65 x 10 14 

100% 

360 ppm 

7.99 x 10 12 

100% 

1.3 x 10 18 

M4 I 

lOOOppm 

4.15 x 10 14 

73.5% 

670ppm 

7.99 x 10 12 

100% 

1.0 x 10 18 

M5 I 

lOOOppm 

1.23 x 10 14 

21.8% 

15000ppm 

7.99 x 10 12 

100% 

4.5 x 10 17 

M6 I 

lOOOppm 

2.42 x 10 13 

4.3% 

15000ppm 

4.04 x 10 11 

5.1% 

1.4 x 10 17 

M7 I 

lOOOppm 

1.67 x 10 13 

3.0% 

15000ppm 

4.96 x 10 10 

0.6% 

6.1 x 10 16 

M8 I 

lOOOppm 

1.34 x 10 13 

2.4% 

15000ppm 

5.72 x 10 9 

0.07% 

2.6 x 10 16 

M9 I 

lOOOppm 

5.87 x 10 12 

1.0% 

15000ppm 

3.17 x 10 9 

0.04% 

1.1 x 10 16 

MUSCLES 

GJ 832 

550ppm 

5.65 x 10 14 

100% 

15ppm 

7.99 x 10 12 

100% 

1.6 x 10 18 

GJ 667C 

330ppm 

5.65 x 10 14 

100% 

0.7ppm 

7.99 x 10 12 

100% 

7.0 x 10 18 

GJ 1214 

1600ppm 

5.65 x 10 14 

100% 

l.lppm 

7.99 x 10 12 

100% 

2.3 x 10 18 

GJ 436 

1600ppm 

5.65 x 10 14 

100% 

4.3ppm 

7.99 x 10 12 

100% 

1.1 x 10 18 

GJ 581 

1400ppm 

5.65 x 10 14 

100% 

2.9ppm 

7.99 x 10 12 

100% 

1.2 x 10 18 

GJ 876 

3400ppm 

5.65 x 10 14 

100% 

6.9ppm 

7.99 x 10 12 

100% 

8.8 x 10 17 


Table 4 

Integrated UV fluxes at the top of the atmosphere (TOA). 


Stellar 

Type 

Ly-a: TOA 

1210-1222A 

ergs cm -2 s -1 

FUV minus Ly-o; 
TOA 1222-2000A 
ergs cm -2 s —1 

FUV TOA 

1000-2000A 

ergs cm -2 s —1 

NUV TOA 
2000-3200A 
ergs cm -2 s -1 

Full UV TOA 
1000-3200A 
ergs cm -2 s _1 

Active 

MO A 

6.9 x 10 2 

2.4 x 10 2 

9.4 X 10 2 

1.4 X 10 3 

2.3 X 10 3 

Ml A 

5.6 x 10 2 

2.0 x 10 2 

7.6 x 10 2 

1.1 X 10 3 

1.9 X 10 3 

M2 A 

4.4 x 10 2 

1.6 x 10 2 

6.1 x 10 2 

8.5 x 10 2 

1.5 x 10 3 

M3 A 

3.5 x 10 2 

1.2 x 10 2 

4.7 x 10 2 

6.5 x 10 2 

1.1 x 10 3 

M4 A 

3.1 x 10 2 

1.1 x 10 2 

4.2 x 10 2 

5.6 x 10 2 

9.8 x 10 2 

M5 A 

1.8 x 10 2 

6.3 x 10 1 

2.5 x 10 2 

3.3 x 10 2 

5.8 x 10 2 

M6 A 

6.4 x 10 1 

2.3 x 10 1 

8.8 x 10 1 

1.2 x 10 2 

2.1 x 10 2 

M7 A 

2.6 x 10 1 

9.2 x 10° 

3.6 x 10 1 

4.8 x 10 2 

8.4 x 10 2 

M8 A 

2.6 x 10 1 

9.0 x 10° 

3.5 x 10 1 

4.7 x 10 2 

8.2 x 10 2 

M9 A 

2.2 x 10 1 

7.6 x 10° 

3.0 x 10 1 

3.9 x 10 2 

6.9 x 10 2 

Inactive 

MO I 

2.7 x 10" 11 

9.4 x 10 -3 

9.5 x 10“ 3 

3.0 x 10 2 

3.0 x 10 2 

Ml I 

5.8 x 10" 12 

2.9 x 10“ 3 

2.9 x 10~ 3 

2.6 x 10 2 

2.6 x 10 2 

M2 I 

1.5 x 10 -12 

9.1 x 10 -4 

9.1 x 10~ 4 

1.7 x 10 2 

1.7 x 10 2 

M3 I 

3.5 x 10 -14 

6.9 x 10 -5 

6.9 x 10 -5 

8.9 x 10 1 

8.9 x 10 1 

M4 I 

6.1 x 10 -15 

2.2 x 10 -5 

2.2 x 10 -5 

5.9 x 10 1 

5.9 x 10 1 

M5 I 

2.9 x 10“ 19 

1.5 x 10“ 8 

1.5 x 10~ 8 

1.3 x 10 1 

1.3 x 10 1 

M6 I 

4.2 x 10" 21 

2.6 x 10“ 9 

2.6 x 10 -9 

3.0 x 10° 

3.0 x 10° 

M7 I 

1.1 x 10" 21 

2.6 x 10~ 10 

2.6 x 10 _1 ° 

1.3 x 10° 

1.3 x 10° 

M8 I 

8.3 x 10" 24 

1.7 x 10~ 10 

1.7 x 10 _1 ° 

5.0 x 10 _1 

5.0 x 10' 1 

M9 I 

1.8 x 10" 24 

1.6 x 10~ 10 

1.6 x 10 _1 ° 

1.3 x 10 _1 

1.3 x 10” 1 

MUSCLES 

GJ 832 

1.5 x 10 2 

3.0 x 10° 

1.5 x 10 2 

3.5 x 10 2 

5.0 x 10 2 

GJ 667C 

1.2 x 10 2 

2.4 x 10° 

1.3 x 10 2 

2.3 x 10 2 

3.6 x 10 2 

GJ 1214 

1.4 x 10° 

9.9 x 10 1 

1.1 x 10 2 

1.8 x 10 2 

2.8 x 10 2 

GJ 436 

6.1 x 10 1 

7.1 x 10° 

6.8 x 10 1 

2.3 x 10 2 

2.9 x 10 2 

GJ 581 

4.2 x 10 1 

6.8 x 10° 

5.0 x 10 1 

2.7 x 10 2 

3.2 x 10 2 

GJ 876 

3.4 x 10 1 

1.1 x 10 1 

4.9 x 10 1 

1.1 x 10 2 

1.6 x 10 2 


N 2 O to a fixed mixing ratio corresponding to the last 


stable run for the inactive model M5. N 2 O is also an 
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Figure 4. Planetary temperature vs. altitude profiles and mixing ratio profiles profiles for H2O, O3, CH4 and N2O (left to right) for an 
Earth-like planet orbiting the grid of active stellar models (top), inactive stellar models (middle), and MUSCLES stars (bottom). Earth-Sun 
profiles are shown for comparison (dashed black lines). 


indirect sink for stratospheric 03 , since about 1 % is con¬ 
verted to NO. Therefore, increasing N 2 0 decreases O 3 
abundance. 

OH concentrations decrease with decreasing UV levels 
(in cooler and inactive M dwarfs) and OH is a primary 
sink for many species in the atmosphere including, but 
not limited to, many biologically interesting species in¬ 
cluding CH 4 , CH 3 , HC1, H 2 , H 2 S and CH 3 C1. CH 3 C1 
concentrations increase for cooler and inactive M dwarfs 
due to decreased stellar UV flux. 

0 2 and C0 2 concentrations remain constant and well 
mixed for all stellar types. 

The planetary surface temperatures range between 297 
- 304K and Bond albedos range between 0.108 - 0.06 for 
M0 to M9 stars, respectively (see Fig. [ 5 b. These albedos 
are lower than Earth’s Bond albedo of (L3 around the Sun 
because the stellar spectral energy distributions peak at 
longer wavelengths for cooler stars where Rayleigh scat¬ 
tering is less efficient, assuming the same total insolation. 

3.1. Effect of Ly-a 

Ly-a is the brightest line in the UV spectra for cool 
stars and accounts for a significant portion of the overall 
UV flux. However, the intrinsic Ly-a flux must be recon¬ 
structed to ac count for interstellar absorption by ne utral 
hydrogen (see Wood et al.|2005 Linsky et al.|2013 ). For 
the observed IvlUyCLES stars, the Ly-a line flux is 13- 
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Figure 5. Calculated Bond albedos for Earth-like planets orbit¬ 
ing the M0-M9 Active (black diamond) and Inactive (red triangle) 
stellar models. 

33% the total UV flux excluding Ly-a. GJ 1214 has no 
observed Ly-a and only an up per limit which is 0.5% 
of the total UV flux (France et al. 2013). In our active 


models, Ly-a ranges from 2 to 17% of the total UV flux 
based on observations of AD Leo. In the inactive models 
(i.e. no chromosphere) Ly-a is negligible compared with 
total UV flux. 
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Temp (K) H 2 0 Mixing Ratio 0 3 Mixing Ratio CH 4 Mixing Ratio N 2 0 Mixing Ratio 

Figure 6. Temperature, H2O, O3, CH4, and N2O mixing ratio profiles where the Ly-a has been increased by a factor of 10 3 , 10®, 10 9 , 
10 42 and 10 15 above the inactive MO model star corresponding to flux levels of 2.7 X 10“ 8 , 2.7 X 10“ 5 , 2.7 X 10“ 2 , 2.7 X 10 1 , 2.7 X 10 4 ergs 
cm' 2 s' 1 , respectively (top) and for GJ 1214 where the Ly-a upper limit (2.4 X 10 —15 ergs cm -2 s -1 ) has been multiplied by a factor of 
10“ 2 , 10~ 1 , 10 1 , 10 2 and 10 3 (bottom) 


We ran two sensitivity tests of the photochemistry of 
Earth-like planets to the amount of Ly-a flux from their 
host star. First we increased the Ly-a by a factor of 
10 3 , 10 6 , 10 9 , 10 12 and 10 15 above MO inactive model 
corresponding to flux levels of 2.7 x 10“ 8 , 2.7 x 10“ 5 , 
2.7 x 10“ 2 , 2.7 x 10 1 , 2.7 x 10 4 ergs cm -2 s _1 respec¬ 
tively. The highest Ly-a value considered is 84x higher 
than that of the MO active model and 180x higher than 
the highest observed Ly-a flux in the MUSCLES stellar 
sample. As seen in Fig. [o] (top), increasing the Ly-a 
flux has only a small effect on the photochemistry un¬ 
til the most extreme case considered where it primarily 
photolyzes N 2 0 and CH 4 . We run our models to around 
60km corresponding to a pressure of 1 x 10 “ 4 bar. The ef¬ 
fect of Ly-a will be mo re p ronounced at pressures lower 
than this as Miguel et al. (2015) has shown for mini- 
neptune atmospheres. 

GJ 1214 is the only MUSCLES star to not have a di¬ 
rectly detected Ly-a flux. An upper limit was placed at 
2.4 x 10 “ 15 ergs cm ' 2 s _1 . We artificially set the Ly-a 
flux to be 10 “ 2 , 10 “ 1 , 10 1 , 10 2 , and 10 3 times this up¬ 
per limit. We tested values both above and below the 
upper limit to see how sensitive an Earth-like planet at¬ 
mosphere is to changes in Ly-a flux when the rest of the 
NUV radiation field is still much higher than the inactive 
model in the previous sensitivity test. As shown in Fig. 
[ 6 ] (bottom), the largest changes are seen in the CH 4 and 
to a smaller extent in stratospheric H 2 O concentrations. 
O 3 and N 2 0 remain relatively constant through the 5 
orders of magnitude change in Ly-a flux. 

Figures [4] and [ 6 ] show it is important to characterize the 
entire UVspectrum including the near UV and the base 
level flux between emission lines, to understand future 
observations of extrasolar planet atmospheres. Ly-a is 
one of the most important lines to characterize. 


4. RESULTS: SPECTRA OF EARTH-LIKE PLANETS 
ORBITING M DWARFS 

Spectra of Earth-like planets orbiting M dwarfs with 
varying UV activity levels show measureable differences 
in spectral feature depths. In the VIS, the depth of 
absorption features is primarily sensitive to the abun¬ 
dance of the species, while in the IR, both the abun¬ 
dance and the temperature difference between the emit¬ 
ting/absorbing layer and the continuum influence the 
depth of features. 

We assume full phase (secondary eclipse) for all spec¬ 
tra presented to sh ow the maximum flux that can be 
observed. Figs 7]|l2 show an Earth-size planet to deter¬ 
mine the specific Hux and planet-to-star contrast ratio. 
A Super-Earth with up to twice Earth’s radius will pro¬ 
vide 4 times more flux and a better contrast ratio. No 
noise has been added to these model planetary spectra 
to provide inputs for a wide variety of instrument sim¬ 
ulators for both secondary eclipse and direct detection 
missions. 


4.1. Earth-like Visible/Near-infrared Spectra (0.4prn - 

4l-im) 

Fig. [7] shows the reflected visible and near-infrared 
spectra from 0.4 to 2/rm of Earth-like planets around the 
M dwarf grid of active, inactive, and MUSCLES stars, us¬ 
ing the SED of the host star. We assume Earth-analogue 
cloud cover. The high-resolution spectra calculated with 
0.1 cm ' 1 steps have been smoothed to a resolving power 
of 800 using a triangular smoothing kernel to show the in¬ 
dividual features more clearly. The Earth-Sun spectrum 
is shown for comparison as a dashed black line. 

Due to the increased stellar flux at shorter wavelengths 
for a G-type star, Rayleigh scattering is much more pro¬ 
nounced for FGK stars, which greatly increases the flux 
from 0.4 to 0.8/rm for an Earth-like planet around a hot¬ 
ter star. Since M dwarfs have stronger NIR emission, and 
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Figure 7. Smoothed, disk-integrated VIS/NIR spectra at the TOA for an Earth-like planet for the grid of active stellar models (top), 
inactive stellar models (middle), and MUSCLES stars (bottom) assuming 60% Earth-analogue cloud coverage model (region 2-4//ni has 
low integrated flux levels and therefore is not shown here). The Earth-Sun spectrum is shown for comparison as a dashed black line. 

at 0.6pnr (the Chappuis band), 0 2 at 0.76pm, H 2 0 at 
0.95pm, and CH 4 at 0.6, 0.7, 0.8, 0.9, 1.0 and 1.7pm. 
Note that shallow spectral features like the visible O 3 
feature would require a very high signal-to-noise ratio 
(SNR) to be detected. 

Fig. [ 8 ]shows details for one of the most notable feature 
in the VIS, the 0 2 A band at 0.76pm in both the rela¬ 
tive flux as planet-to-star contrast ratio (left), and the 
reflected emergent flux for a 60% cloud cover model for 
M0-M9 model stars (top middle), and for the six MUS¬ 
CLES stars (top right). Note that the detectability of 
the 0 2 feature in reflected light is similar for active and 
inactive models since the stellar flux at 0.76pm is activ¬ 
ity independent. The relative flux shows an equivalently 
deep feature for each case due to a constant mixing ra¬ 
tio of 21%. However, the oxygen feature in absolute flux 
(Fig. 8 , upper middle panel) becomes increasingly diffi¬ 
cult to detect for the later stellar types. For the latest 
M stellar types modeled here, the detection of the 0 2 
feature requires a very high SNR to detect even if the 
planet has an active photosynthetic biosphere like the 
Earth (see Fig. [f|. If one assumes a blackbody radia¬ 
tion Planck function, rather than a realistic stellar model, 
when calculating the shape of the reflected light curve, 
the reduction in the feature’s depth on moving to later 
stellar types (i.e. going from blue to red lines in Fig. [ 8 j 
2nd and 3rd columns) is less pronounced. The 0 2 fea- 



Wavelength (jim) 


Wavelength (|xm) 


Wavelength (nm) 


Figure 8. O 2 feature at 0.76jLtm relative flux as planet-to-star 
contrast ratio (left), and the reflected emergent flux for a 60% 
cloud cover model for M0-M9 model stars (top middle), and for 
MUSCLES stars (top right). Bottom row shows the same but with 
a blackbody used to calculate the absolute flux. The observable 
O 2 feature becomes increasingly difficult to detect for later stellar 
types due to the decreasin g s tellar flux in this wavelength region. 
Coloring is same as in Fig. [7] 

the 1-2/mi flux is larger around Earth-like planets orbit¬ 
ing M dwarfs than for the Earth-Sun equivalent. The 
most notable features in the VIS/NIR spectra are O 3 









































UV Radiation on Earth-like Planets Orbiting M dwarfs 


11 


ture is pronounced for all MUSCLES stars since none of 
those stars have a low enough T e g for the O 2 feature to 
be diminished by the spectral energy distribution (SED). 
GJ 876 is the coolest MUSCLES star with T e g= 3129K 
and in our models the O 2 feature becomes most obscured 
for stars with T e g = 2300K - 2600K. 

CH 4 also has several features of interest in the 
VIS/NIR range at 0.6, 0.7, 0.8, 0.9, 1.0 and 1.7pm (see 
Fig. |9b. The CH 4 feature is deeper for less active and 
coolerM dwarfs. In particular, the CH 4 features at 0.8, 
0.9 and 1.0 pm become much more pronounced for the 
cooler M dwarfs and for some of the MUSCLES stars 
such as GJ 876. 



1.5 1.6 1.7 1.8 1.9 1.5 1.6 1.7 1.8 1.5 1.6 1.7 1.8 

Wavelength (urn) Wavelength (|im) Wavelength (|xm) 


Figure 9. CH 4 features at 0.6, 0.7, 0.8, 0.9, and 1.0/mi (top row) 
and the CH 4 feature at 1.7/mi (bottom row) for the relative flux as 
planet-to-star contrast ratio (left), and the reflected emergent flux 
for a 60% cloud cover model for M0-M9 model stars (middle), and 
for MUSCLES stars (right). The observable CH 4 feature becomes 
deeper for planets orbiting cooler and less active star models. Col¬ 
oring is same as in Fig. m 


H 2 O several features in the VIS/NIR at 0.95, 1.14 ,1.41 
and 1 . 86 pm and are deeper for cooler and less active stars 
(see Fig. [ 7 ]). 

Clouds increase reflectivity but can also decrease the 
depth of all observable features compared to the clear sky 
case. It will be difficult, therefore, to remotely determine 
the absolute abundance of a molecule without a well- 
characterized temperature vs. pressure distribution as 
well as cloud profile. 


4.2. Earth-like Infrared Spectra, IR (4pm - 20pm) 


Fig. 10 shows the thermal emission spectra from 


4 to 20pm of Earth-like planets with Earth-analogue 
cloud cover around the M dwarf grid of active, inac¬ 
tive, and MUSCLES stars using the stellar SED. The 
high-resolution spectra have been smoothed to a resolv¬ 
ing power of 150 using a triangular smoothing kernel to 
show the resolution expected by JWST. The Earth-Sun 
IR spectr um is shown for comparison as a dashed black 
line. Fig. [lT] shows the individual component gas contri¬ 
butions ofthe dominate gases (H 2 0, CH 4 , C0 2 , CH 3 C1, 
0 3 , and N 2 0) to the final IR planetary spectrum for an 
M9 active and inactive stellar model. 

The depth of the O 3 feature at 9.6pm decreases for 
planets orbiting cooler and less active M dwarfs, as ex¬ 


pected due to lower O 3 abundances for lower UV incident 
flux and also due to the decreased temperature differ¬ 
ence between the O 3 emitting layer (around 40km) and 
the surface temperature, which is larger for the earlier M 
dwarfs as seen in Fig. [4j 

The CH 4 feature at 7.7pm, decreases in depth for 
planets orbiting cooler M dwarfs despite increasing CH 4 
abundances making it difficult to remotely determine the 
CH 4 abundance without a well-characterized tempera¬ 
ture vs. pressure profile. Note that the 7.7pm feature 
is partially obscured by the wings of the H 2 0 feature at 
5-8pm. 

The N 2 0 features are the most striking addition in 
the inactive models with low UV flux. N 2 0 has features 
at 7.75pm (overlapping with the CH 4 feature), 8 .5pm 
10.65pm and 16.9pm (see Fig. 
with decreasing UV flux. In t 


11 ) that become deeper 
le active M0-M9 stellar 


models we do not see any strong N 2 0 features, although 
it does contribute to the depth of the 7.7 CH 4 feature. 
For the inactive M5-M9 models we see a strong N 2 0 fea¬ 
ture at 10.65pm and wide absorption from 14-19pm due 
to the strong build-up of N 2 0 in low UV environments as¬ 
suming an Eartli-like biological surface flux of 7.99 x 10 12 
g yr ^ 1 (see §2 and §3). However, even for modest mixing 
ratios of 15ppm to 34ppm in GJ 832 and the M0 inac¬ 
tive models, respectively, we see a strong N 2 0 feature at 
16.9pm. These concentrations are not much higher than 
the model mixing ratios calculated for the other MUS¬ 
CLES stars and the later active M dwarf models. N 2 0 
has been considered to be a strong bio signature (see e.g. 
Segura et al.||2005~ Seager et al.||2012 ) and will be easier 


to detect around M dwarfs - especially inactive ones - 
than around FGK stars. 

The C0 2 absorption feature at 15pm does not have a 
central emission peak, usually seen in F and G stars, due 
to the more isothermal stratospheres of planets around 
all M dwarfs. 

The depth of the H 2 0 feature at 5-8pm decreases for 
later M dwarfs despite increasing stratospheric H 2 0 con¬ 
centrations. 

Clouds reduce the continuum level and the depth of 
all of the observable features compared to the clear sky 
case. Therefore, it will be difficult to remotely deter¬ 
mine absolute abundances without a well-characterized 
temperature vs. pressure distribution as well as cloud 
profile. 

Observability of Biosignatures: Detecting the combi¬ 
nation of 0 2 or O 3 and a reducing gas like CH 4 in emer¬ 
gent spectra and in secondary eclipse measurements re¬ 
quires observations either in the IR between 7 and 10pm 
that includes the 7.7pm CH 4 and 9.6pm O 3 features, 
or observations in the VIS to NIR between 0.7 to 3pm 
that includes the 0.76pm 0 2 and 2.4pm CH 4 features in 
the observed spectral range. The strengths of absorption 
features (see Fig. [ 7 ] and fll| for VIS/NIR and IR, respec¬ 
tively) depend on the effective temperature of the host 
star as well as its UV flux and vary significantly between 
stellar types (see Fig. HD.- . 

In the IR, CH 4 at 7. 7pm is more easily detected at low 
resolution for earlier grid stars (M0-M4) than for later 
grid stars (M5-M9) in the active models. 

The 9.6pm O 3 feature is deepest in earlier and/or more 
active M dwarfs where there are more UV photons. For 
M5-M9 inactive models, O 3 is not detectable in low res- 
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Figure 10. Smoothed, disk-integrated IR spectra at the TOA emitted by an Earth-like planet orbiting the grid of active stellar models 
(top), inactive stellar models (middle), and MUSCLES stars (bottom) assuming 60% Earth-analogue cloud coverage model. The Earth-Sun 
spectrum is shown for comparison as a dashed black line. 


olution at 9.6/im due to the extremely low UV flux for 
these models. The narrow O 2 feature in the VIS at 
0.72/xm becomes weaker for cooler M dwarfs even though 
the mixing ratio of O 2 remains fixed in our simulations at 
21%. This is a consequence of the faint spectral energy 
distributions at shorter wavelengths for both active and 
inactive late M stellar models. 

For the modern Earth, N 2 O and CH 3 CI do not con¬ 
tribute substantially to the spectrum due to their low 
mixing ratios and will likely be undetectable by the first 
low-resolution and photon-limite d exo pl anet atmosp here 
chara cterization missions ( Selsis|2000| Kaltenegger et al. 
2007). However, both N 2 (J and CH 3 C1 reach detectable 
levels in the IR spectra for our models of cooler and less 
active stars due to low photolysis rates even thou gh N 2 Q 
has many absorption features in the IR (see Fig. 11 1 . 
N 2 0 is considered a strong bio signature b ecauseu iere 


are no significant abiotic sources (Des Marais et al. 2002). 
One MUSCLES star, GJ 832, has a noticeable 1N 2 U tea- 
ture at 17/im with 15ppm for a modern Earth-like flux. 
This case is interesting because GJ 832 has the highest 
total FUV flux at wavelengths where N 2 O is photolyzed 
(1000-2400A). We find that the high Ly-a flux of GJ 832 
promotes the destruction of O 3 , producing more 0( 1 D) 
which then reacts with N 2 to form N 2 0. 

For planets orbiting inactive M5-M9 star models we 
observe a sharp increase in N 2 O concentrations as a con¬ 


sequence of few available UV photons. N 2 0 features 
dominate the spectrum for these late, inactive M dwarfs 
and can be detect ed in the IR at 4-5/rm, 8 - 11 /xm, and 
16-19/xm (see Fig. 111. In an atmosphere dominated by 
N 2 O, such strong absorption throughout the IR could 
obscure the 7.7//m CH 4 feature. It is unknown whether 
a strong build-up of biotic N 2 0 would be physically pos¬ 
sible around stars with little or no chromospheric flux. 
Given the existence of a few GA LEX stars w ith l ittle ex- 
cess chromospheric NUV flux (Shkolnik & Barman |2014 1 
and only an upper limit established for the Ly-a flux from 
GJ 1214, a small fraction of M dwarfs could exhibit low 
UV fluxes. However, the MUSCLES data set shows that 
all six observed M dwarfs have sufficient UV flux to pho- 
tolyze N 2 O and to prevent “run-away” N 2 O build-up. 

CH3CI contributes to our IR spectrum from 13-14yLtm 
in the short wavelength wing of the CO 2 and N 2 O fea¬ 
tures and could also be detectable, depending on the CO 2 
and N 2 0 concentrations present (see Fig. [Til. 

For clear sky models, the vegetation reaedge (VRE) 
surface feature is detectable in low resolution spectra 
due to the order of magnitude increased VRE reflectance 
between 0.7/im and 0.75/xm for all M grid stars assum¬ 
ing that the e xoplanets of these host s tars have similar 


plant life (see Rugheimer et al. 2013 for FGK stars). 


Clouds partly obscure this feature compared to the clear 
sky case, although the increase in flux can be seen at 
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Figure 11. Individual spectral components for H 2 O, CH 4 , CO 2 , CH 3 CI, O 3 , and N 2 O comprising the full IR spectra (shown by a dashed 
black line) for an Earth-like planet orbiting a M9 active stellar model (top) and M9 inactive stellar model (bottom). 


0.7/irn in Fig. [ 7 ] which include Earth-like clouds. Due to 
the shift in available photons to longer wavelengths for 
M dwarfs, a different photosynthesis biochemistry could 
have evolved, resulting in a di fferent b ut potentially ob¬ 
servable vegetation signature ( Kiang et al.| [2007). 

Observability of Spectral Features Note that we have 
not added noise to these model spectra in order to be use¬ 
ful as input models for a wide variety of instrument sim¬ 
ulators for both secondary eclipse and direct detection 
simulatio ns. Different instrument simulators for JWST 
(see e.g. |Deming et al.|2009; Kaltenegger & Traub|2009 1 
explore the capability of J WST’s MIRI and JNIRspec in- 
struments to characterize extrasolar Earth-like planets 
for near-by as well as luminous host stars. Several groups 
are providing realistic instrument simulators that can 
be used to determine the detectability of these absorp¬ 
tion features. Future ground and space based telescopes 
are being designed to characterize exoplanets as small as 
Earth-like planets and will provide opportunities to ob¬ 
serve atmospheric features, especially for Super-Earths 
with radii up to twice Earth’s radius and therefore four 
times the flux and planet-to-star contrast ratio levels for 
Earth-size planets as shown in Fig. [12] 

In addition to measuring the size ofthe planet, future 
observations will occur at different phases throughout 
the planet’s orbit. The maximum observable planetary 
flux in the visible spectrum scales with the illuminated 
fraction of the planet that is visible to the observer. At 
quadrature, representing an average viewing geometry, 
the contrast ratios presented in Fig. [12] will be a factor 
of ~ 2 lower in the visible. In the IR, the maximum flux 
remains constant throughout the planet’s orbit, assum¬ 


ing a similar temperature on the planet’s day and night 
side. 


5. DISCUSSION 

Our active and inactive stellar models represent the 
extreme ends of stellar activity for M dwarfs. AD Leo, 
our active star proxy with H a in emission, was chosen to 
represent a young, active star. Our inactive stellar mod¬ 
els assume no chromospheric activity, which is the lowest 
level possible of UV flux. More observations are needed 
to determine the lowest possible level of UV flux around 
an M dwarf, which still might be significantly higher than 
the photosphere only models. Many M dwarfs exhibit 
a wide range of activity. Therefore, we use the MUS¬ 
CLES star database to represent observed M dwarfs with 
a range of activity between the two extremes. The sim¬ 
ulated planets orbiting the MUSCLES stars share more 
similarities with the planets orbiting the active stellar 
models particularly in terms of O 3 , CH 4 , and N 2 O con¬ 
centrations, which is to be expected since they have UV 
fluxes detectable with HST. 

Ly-a is an important line to characterize because it is 
the strongest line in the FUV part of an M dwarf spec¬ 
trum. We note that in addition to properly reconstruct¬ 
ing the Ly-a line, determining the base level of flux and 
other emission lines in the FUV and NUV are also impor¬ 
tant. This is because Ly-a by itself does not significantly 
impact the photochemistry of an Earth-like planet atmo¬ 
sphere below about 60km unless the Ly-a flux exceeds 
the highest observed levels observed for the MUSCLES 
stars (see §3 and Fig. [ 6 j> . 

Our present models assume a steady state. Young M 
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Figure 12. Contrast ratios for Earth-like planets orbiting M0-M9 stars in full illumination. Active model contrast ratios are shown here. 
The contrast ratios of spectral features of planets orbiting inactive M dwarfs have similar levels as active M dwarfs in our models. The 
black dotted lines show contrast ratios for the Earth orbiting the Sun. 


dwarfs exhibit strong flares which may impact the at¬ 
mosphere on timescales relevant to the photochemistry. 
Future work will consider activity ind uced variabilit y in 
an exoplanet’s atmosphere. See Segura et al. (2010) for 
models on time-dependent behavior of biosignafures for 
flaring M dwarfs. 


6. CONCLUSIONS 

We show spectral features for terrestrial atmosphere 
models of Earth-like planets from the VIS to the IR for 
planets orbiting a grid of MO to M9 host stars (T e g = 
2300K to 3800K) using active stellar models, inactive 
stellar models, and the sample of M dwarfs with observed 
UV fluxes. Our grid comprehensively covers the full M 
stellar range and activity from a proxy of an active flare 
star, AD Leo, to the lower limit of activity with no chro¬ 
mospheric contribution to the UV. 

We discuss the atmospheric model results in §3 and 
the detectable spectral features with a focus on how the 
UV flux environment affects the emergent spectra in §4. 
Increased UV flux environments in M dwarfs is primar¬ 
ily a property of younger and earlier type stars. Higher 
UV environments produce: increased concentration of 
O 3 from photolysis of O 2 , increased O( 1 D) from pho¬ 
tolysis of O 3 , increased OH from reaction of 0( 1 D) and 
H 2 0, decreased stratospheric H 2 0, and decreased CH 4 , 
CH 3 CI, and N 2 0 from photolysis and reactions with OH 
(see Fig. 41). Very inactive M dwarfs may have extremely 
low UV flux levels as suggested by GALEX observations 
in the NUV. More observations are needed to determine 
the lower limit of fluxes in the FUV, in particular the 
Ly-a emission. Such low levels of UV flux could lead to 
a build-up molecules such as CH 4 and N 2 0 in the atmo¬ 
sphere of an Earth-like planet, with N 2 0 being a strong 
biosignature due to a lack of photochemical sources. 

Terrestrial planets orbiting early M spectral type stars 
are the best targets for observing biosignatures, such as 
0 2 or O 3 in combination with a reducing species like 
CH 4 . Note that the 0 2 line becomes increasingly dif¬ 
ficult to detect for late M dwarfs (see Fig. [ 8 ]). While 
not observable for an exoplanet around an FGK star, 
N 2 0 may be observable at 10.7/rm or at 17/xm at Earth- 
level emissions for planets orbiting M dwarfs (see inactive 
models and model for GJ 832). 


Our results provide a grid of atmospheric compositions 
as well as model spectra from the VIS to the IR for JWST 
and other future space and ground-based direct imaging 
and secondary eclipse missions in terms of instrument 
design and will help to optimize their observation strat¬ 
egy. The model spectra in this paper are available at 
www.cfa.harvard.edu/~srugheimer/Mspectra/. 

ACKNOWLEDGMENTS 

We especially thank Andrew West and Evgenya L. 
Shkolnik for useful conversations about M dwarfs. We 
also would like to thank Kevin France for discussions 
concerning the MUSCLES database. This work has 
made use of the MUSCLES M dwarf UV radiation field 
database. We would also like to acknowledge support 
from DFG funding ENP KA 3142/1-1 and the Simons 
Foundation (290357, Kaltenegger). 

The UV radiation field for AD Leo was obtained from 
the Multimission Archive at the Space Telescope Science 
Institute (MAST). STScI is operated by the Association 
of Universities for Research in Astronomy, Inc., under 
NASA contract NAS5-26555. Support for MAST for 
non-HST data is provided by the NASA Office of Space 
Science via grant NAG5-7584 and by other grants and 
contracts. 


APPENDIX 


GJ 667C: GJ 667C is an M3-4 dwarf at 6.8 pc (van 
Leeuwcn 2007) with an age esti mate greater than 2 Gyr 
(|Anglada-Escude et al.| |2013b|). GJ 667C has a T c g 
= 3330K. Fe/Hj OT55, and M = O.33OM 0 ± 0.019 
(Anglada-Escude et al. 2013b). Previous estimates had 
the effective temperature 350K higher and the spectral 
design ation to be Ml 5 due to assuming a higher metallic - 
ity (see|Geballe et al.|2002||Anglada-Escude et al.|2013b|). 
We use a radius of R = 0.35/b., consistent with the lumi - 
nosity, L = 0.01370 L 0 ((Anglada-Escude & Tuomi 2012|). 
We merged a PHOENIX UT-Settl spectrum for l e ff = 
3350K, [Fe/H] = -0.55, log(g) = 5, and v sin i = 3 km 
s -1 with the MUSCLES UV spectrum at 2800 A. GJ 
667C has less wavelength coverage than the other MUS¬ 
CLES stars. We approximated the rest of UV radiation 
field from the FUV and NUV spectral energy distribu¬ 
tion of GJ 832, which has a similar spectral type (France 
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et al. 2013), scaled to the distance and Ly-a and Mg II 
emission line strength of GJ 667C. 


GJ 832: GJ 832 is an M dwarf at 4.95 pc (van Leeuwen 


2007) with no age estimate and is the least well charac¬ 


terized star in the MUSCLES program sample. It has 
a T a a. =_3620K (NStED value interpol ated as des cribed 
Bessell||1995 ) , R = O.48R 0 ( Joh nson fe Wright||1983 ), 


[Fe/H| = -0.12 ( Johnson fe Apps|[2009 ), and log g = 4.7 


(|Bailey et al. 2008|). We merged a PHOENIX BT-Settl 
spectrum with 1 e ff = 3620K, [Fe/H] = -0.12, log g = 4.7, 
and a v sin i = 3 km s -1 with the HST UV spectra from 
the MUSLES program at 2800 A. 

GJ 1214: GJ 1 214 is an M6 dwarf at 14.55± 0.13pc 
(Anglada- Escude et al. 2013a ) with a n age estimate of 
6+3 Gyr (Charbonneau et al. 2009). A more accu¬ 
rate parallax measurement increased the previous dis¬ 
tance estimate by 10% and thus the luminosity and mass 


has been shifted from previous values as well (Anglada- 


Escude et al. 2013aI. With the new parallax and lumi¬ 
nosity, GJ 12 14 has a radius of 0.211 Rr. : . T R ff = 3250K, M 
= 0.17 6 Mp) ( Anglada-EscudbetalJgOlSa ), and [Fe/H] = 


+0.05 ( Anglada-Escude et al.||2Qj3a ~ Neves et al.||2012|). 


=T" 


et al. 2010; Delfosse et al. 1998; Reiners & Basri 20081 

West & Basri 

2009) and log g = 4.991 (Charbonneau 

et al. 

2009). 

We merged a PHOENIX BT-Settl spec- 


and v sin i — 1 km s _1 with the MUSCLES UV spectrum 
at 2800 A. GJ 1214 is the only MUSCLES star to have 
a non-detection of Ly-a, and t hus an up per limit is used 
and discussed in depth in France et al. (2013). 


GJ 436: _ GJ 436 is an M3 dwarf at 10.1 pc (van 


Leeuwen 2007b with an age estimate of 6.5-9.9 Gyr 


Satie et al. 2005). GJ 436 has a T eff = 3416K, R = 
'.455^~0UT8 A 0 and M = 0.507 (+0.071/-0.062) M 0 
(|von Brau n et al. 2012). GJ 436 has solar me tallicity, 
|Fe/H|=0 (|Roj as- Ayala et al.|2010 ), log g = 5.0 flManess 


et al. 2007[) and a v sin * <1 km s 1 (Marcy & Chen 
1992p . We merged a PHOENIX BT-Settl spectrum with 


Jeff = 3416K, [Fe/H] = 0. log g = 5 and a v sin i = 1 
km s _1 with the HST UV spectra from the MUSCLES 
program at 2800 A. 

G1 581: G J 581 is a M3 dwarf at 6.2± O.lpc (Ivan 
Leeuw en|2007|) with age estimate of 7-11 Gyr (|Seisis et al. 
2007|. GJ 581 has a T eS = 3498.0K+ 56.0K, R = 0.299 


± 0.0 10 i? Q , and log g = 4.96+0.08 (von Braun et al. 


20111. GJ 581 has a metallicity s lightly subsolar, |Fe/H| 


= -0.02 ( Rojas-Ayala et a L||2010[) and an upper limit on 
v sin * < 2.1 km s _1 (Dclfosse et al.|l998). We merged a 
PHOENIX BT-Settl spectrum with T eS = 3498K [Fe/H] 
= - 0.02, log g = 4.96, and v sin i = 2.1 km s -1 (Correia 


et al. 2010) with the MUSCLES UV spectrum at 2800 
- 


GJ 876: GJ 876 is an M d warf at 4.69pc (van Leeuwen 
20071 with an age 0.1-5Gyr ( Correia et al.|2QTo[ ). GJ 87(7 


has Jfiff = 3129+19 K and R = 0.3761 ± 0.0059 Rq (von 


Braun et al. 2014). We merged a PHOENIX B T-Settl 


spectrum with 2 (A = 3129K, [ Fe/H] = +0.19 ( Rojas- 
Ayala et al.||2012 ), log g = 4.89 ([B ean et aT 12006 ), and 


sin i = 1.38 km s 1 (Correia et al. 2010[ ) with the 
MUSCLES UV spectrum at 2800 A. 
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